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ABSTRACT 

The newly established luminosity functions of high-z galaxies at 4 < z < 10 can provide a stringent 
check on dark matter models that aim to explain the core properties of dwarf galaxies. The cores 
of dwarf spheroidal galaxies are understood to be too large to be accounted for by free streaming of 
warm dark matter without overly suppressing the formation of such galaxies. Here we demonstrate 
with cosmological simulations that wave dark matter, i/)DM, appropriate for light bosons such as 
axions, does not suffer this problem, given a boson mass of > 1.2 x lCV 22 eV (2a). In this 
case, the halo mass function is suppressed below ~ 10 10 M 0 at a level that is consistent with the 
high-z luminosity functions, while simultaneously generating the kpc-scale cores in dwarf galaxies 
arising from the solitonic ground state in ^DM. We demonstrate that the reionization history in 
this scenario is consistent with the Thomson optical depth recently reported by Planck, assuming a 
reasonable ionizing photon production rate. We predict that the luminosity function should turn over 
slowly around an intrinsic UV luminosity of M-yy 16atz>4. We also show that for galaxies 
magnified >10x in the Hubble Frontier Fields, ^>DM predicts an order of magnitude fewer detections 
than cold dark matter at z > 10 down to -Muv ~ —15, allowing us to distinguish between these very 
different interpretations for the observed coldness of dark matter. 

Subject headings: cosmology: theory - dark ages, reionization, first stars - dark matter - galaxies: 
abundances - galaxies: evolution - galaxies: high-redshift 


1. INTRODUCTION 

The transitional stage when density perturbations first 
collapse to become galaxies is now being reached with 
data of unprecedented depth. New polarization measure¬ 
ments with the Planck satellite indicate an optical depth 
of only r e = 0.066 ± 0.016, for electron scattering of the 
cosmic microwave background (CMB), which translates 
into an approx imate i n stantaneous rei onization redshift 
of z ~ 8.8l 0 4 (Planck Collaboration et ah 2015). 

This may be compared with the ionization state of hy¬ 
drogen absorption lines revealed by distant, bright quasi- 
stellar objects (QSOs) and gamma-ray bursts (GRBs), 
which has demonstrated that inter-galactic hydrogen has 
remained highly ionized since at least z ~ 6. Be- 


yond which evidence of a patclr 
some QSO and GRB spectra ( 

f opacity is claimed in 

Chornock et al. 

2014; 

Meiandri et al. 

2015 

Becker et 

al. 2015) with a 

small 


mean enhancement in the av erage neutral fraction by 


z = 6.5 (Hartoog et al. 2014). The forest at 2 > 6.5 is 
not yet measured with any bright source, but a large 
jump in neutral hydrogen fraction may be implied in 
the range 6.5 < 2 < 7.0 from the statistical absence of 
strong Lyman-a emission f rom the most distant galaxies 


amenable to spectroscopy (Pentericci et al. 2014 Stark 
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et al. 2015). Exceptions have been found with strong 
Lyman-a emission at z ~ 7.7 (Oesch et al. 20151, and 
also the currently highes t redshift galaxy e stablished by 
spectroscopy at 2 ~ 8.7 (Zitrin et al.|2015). A relatively 
sharp reionization transition is considered likely if galax¬ 
ies are the dominant source of reionization, based on ad¬ 
vanced hybrid techniques (Mesinger et al. |2015 1 which 
justifies the instantaneous re ionizati on r edshift approxi¬ 
mation (Dijkstra|2014 Mitra et al.|2015 ). More detailed 
3D modeling of this transition remains very challenging 
when incorporating all relevant proces ses w hich may af- 


fect ultraviolet (UV) ionization (e.g. , Springel fe Hern - 

F 


quist 12003 Bromm & Yoshida 2011 Wise et aT | 2U14 ) 


including early galactic outflows (eTg., lfye_etaLJ 20U2| 


and their feedback effects (e.g., Scan napicco fe Broad- 
~| |2007j Booth et al. 2012|77 
aiue of T P implies that galaxi 


hurst 2001 Pieri et al. 


‘This relatively low value of r e implies that galaxies may 
not be expected to be found in abundance at redshifts 
much higher than z ~ 9. Indeed beyond this redshift 
only a handful of galaxies are claimed in the deepest 
fields, among which the current most distant galaxy has 
10.7, which is discovered in the CLASH prog ram and 


is highly magnified by gravitational lensing (Coe et al. 


2013). Two other reliably estimated high-z galaxies are 


also known behind the new, highly magni fying Hu bble 
Front ier Field (HF F) clusters at z ~ 9.6 (Zheng et al. 


20121) and z ~ 9.8 (|Zitrin et al. 20141. Nothing beyond 
this has been found in HIT 1 so far , despite the high mag¬ 


nifications by these clusters flZitrin fe Broadhurst 2009 


Lam et ad7||2014 Diego et al.||2015[) and the great depth 


of these images in the near infrared, with their potential 
to detect even higher redshift galaxies to a photometric 
limit of z ~ 11.5 (Coe et al. 20151). 

The UV luminosity function (LF) of distant galaxies 
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is now well constructed at z ~ 4 — 10 relying mainly 
on dropout galaxies de tached in deep field searches (e.g., 


Bouwens et al.||2015b hereafter B15b). The LF is seen 
to steadily evolve to low number densities of high red- 
shift galaxies and to steepen at the faint-end slope as it 
does so (B15b). Evolution is also seen in terms of the 
mean sizes of dropout g alaxies which steadily decrease 
with increasing redshift (Bouwens et al.|[2003; Holwerda 


et al. 20151. Currently, the behavior of the U V lumi¬ 


nosity density at 2 > 8 is hotly debated with evidence 
of an accelerate d decline at 2 > 8 by B15b but with a 


counter claim by McLeod et al. (2015) for the HFF with 
a reliance on parametr ic lens mod eling. This latter claim 


is at odds with Coe et al. (20151 who provide some evi¬ 
dence of a deficit at high redshifts for the HFF, based on 
the first two completed clusters. Of course in this 2 > 9 
redshift range data is restricted to fewer detections in 
only infrared passbands lying close to magnitude limits, 
so that crucial conclusions regarding galaxy formation 
are still uncertain. 

In this paper, we examine the high-redshift galaxy for¬ 
mation in the c ontext of a wave d ark matter model, 


known as ^DM (Schive_et ah 20 14a hereafter SCB14a) 
or fuzzy dark matter {Hu et aT° 20o01. In this scenario, 


dark matter is assumed to be composed of extremely 
light bosons, such as axion-like part icles proposed by 


strin g theory (Arvanitaki et al. 2010) or non-QCD ax- 


ions (Chiueh 2014). They are non-thermally generated 


((Turner 

19831 Goodman 2000 

Bohmer & Harko 

2007 

Sikivie 

Yang||2009; Davidson 

20151 Guth et al. 

2015) 


When self-interaction is negligible, the evolution o 


is governed by the Schrodinger-Poisson equation ( Rufhni 

fc Bonazzola|19 69 ; Seid el fe Suen|1990||Widrow fe Kaiser 


1993 


,Hu et al.pOOOt |Woo fc Chiuehfj2009| |Schive et al. 


2014b), with a single free parameter, to^,, the dark mat- 
Lef particle mass. 

The most prominent feature in ^DM is that the un¬ 
certainty principle counters gravity below a Jeans scale, 
resulting in a suppression of halos below ~ 10 10 M 0 and a 
flat density profile within ~ 0.1 —1.0 kpc of the centers of 


(jKhlopov et al. 1985[ Peebles 2000 Hu et al. 2000 

Matos 

& Arturo Urena-L6pez||2001; Lee & Lim|2010; M 

iarsh & 

Ferreira]2010; Schive et al.|2014a|). This boson mass scale 

can naturally arise in a non-QCD axion model ( 
20141, lending support for the very light bosor 
ADM model has become a viable dark matter ca 

Chiueh 
1 . The 

ididate 

(e.g., Woo & Chiueh 2009 Mielke & Perez 2009 Chava- 


nis 2011; Suarez & Matos 2011:; Robies & Matos 

2012 

Marsh & Silk 2014 Rindier-Daller & STiapiro 2014 

Lora 

& Magana 2014) Bray & Goetz| 2014 

Suarez et al. 

2014 

Bozek et al. 2015; Marsh & Pop 2015 

Martinez-IV 

edina 


et al.||2015MCluzmjn fc Lora-Claviio||26l5| |Madarassy fc 


Toth 2015 [ |Harko fc Lobo||2015[ ), especially given the in¬ 


creasingly strict limits of non-detections of the weakly 

in the standard 


cold dark matter (CDD 
observable properties 0 

d, Akerib et al. 2014). Various 

this model have been proposed 

1 Amendola & Barbieri 

J2006 Arvanitaki & Dubovsky 

2011 

Schive et al. 2D 

l4a||b; Khmelnitsky & Rubakov 

2014 

liozek et al. 12015 

Van Tilburg et al. 20151 IStadnik 

& Flambaum 2015|). 


the )/)DM model have recently generated exciting results 
(SCB14a). We have directly demonstrated that indeed 
the large-scale structure of i/jDM is statistically indis¬ 
tinguishable from CDM, but differs radically on small 
scales, where t/jDM halos form central solitonic cores 
surrounded by fine-scale, large-amplitude granular tex¬ 
tures. By applying a Jeans analysis to the stellar phase- 
space distribution in the Fornax dwarf spheroidal (dSph) 
galaxy, which is know n to have a distinct core (e.g., 


Amorisco et al. 2013), we determine TO 22 = 0.8 ± 0.2 
(ler), thereby providing the crucial normalization of this 
model (SCB14a). 

From our numerical simulations and theoretical argu¬ 
ments based on the scaling symmetry of the Schrodinger- 
Poisson equation and the uncertainty principle, we sub¬ 
sequently _derived_a unique core-halo mass relation in 
V'DM ( |Schive et al. 2014b. hereafter S14b), M c oc (1 + 

2)V 2 M7V where M c and M h are the core mass and 
halo mass, respectively, and 2 is redshift. This relation 
predicts that massive galaxies with M h ~ 10 12 M 0 at 
2 ~ 8 will have compact solitons of M c ~ 10 9 Mq within 
~ 60 pc. Our simulations show that these dense solitonic 
cores form promptly after halo collapse, and t hus m ay 
help to ex plain the early onset of QSO activity (Trakht- 


enbrot et al. 20151 by acting as a massive focus lor gas 


accretion. 

Another key prediction of i/)DM is that galaxy for¬ 
mation is delayed relative to CDM because of the in¬ 
herent Jeans scale. The preliminary results of SCB14a 
showed that the first galaxies form at 2 ~ 13 with 
Mh ~ 10 9 — 10 10 Mq, assuming 77122 ~ 1.0 fixed by the 
scale of clSph galaxy cores as described above. Halos 
below ~ 10 9 Mq are significantly suppressed. We stress 
that the particle mass, 77122 , is the only free parameter 
here assum ing that the dark ma tter is made entirely of 


t/jDM (see Marsh & Silk 
ipDM model). 


2014 for a mixed CDM and 
er the TO 22 , the greater the 
Here we aim to 
1 . 0 ) 


The smal 

difference between CDM and t/iDM. 
establish whether a similar particle mass (77722 
can satisfy both the observed properties of dSph galax¬ 
ies and the constraints from high-z observations, such as 
galaxy counts, reionization history, an d Lyman-q forest , 
although some tension seems to exist (Bozek et al.|2015). 
It is a well known issue for warm dark m atter (WDM), 
usual ly termed as the Catch 22 problem (Maccio et al. 


2012), where the kpc-scale cores in dSph galaxies require 


too small a WDM particl e mass that is in contradiction 
with high-z observations ( Schneider et al.| 2014 Schultz 
et al.|2014||Lovell et al.|2014J). Since the relation between 


core radius and power sp ectrum suppress i on are differ¬ 


ent in ADM and WDM (Hu et al. 


2000 Schive et al. 


2014b Marsh & Pop 20151, in this work we examine m 


detail wlifiHieF^nTRFis’cIear of this serious problem fac¬ 
ing WDM. 

In this paper we conduct cosmological simulations to 
study the evolution of halo mass function (MF) in the 
■0DM scenario, and connect it to the recently established 
galaxy UV LF at 4 < 2 < 10. We explore the results of 
different t/jDM particle masses ranging from 77122 = 0.8 
to 3.2. We predict the evolution of LF beyond the cur¬ 
rent observational limit as a future test to distinguish 
between CDM and t/jDM. We also perform analytic cal¬ 
culations to study the reionization history in this context, 
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and compare it to the Thomson optical depth recently 
reported by Planck. All magnitu des in this paper are 
quoted in the AB system (Mab, Oke fe Gunn||1983 ). 

The paper is structured as follows. In Section [2] we de¬ 
scribe our simulation setup, including initial conditions 
and simulation characteristics. We show the -0DM halo 
mass function in Section [31 and compare it with obser¬ 
vations in Section [1J Finally, we discuss and summarize 
our results in Section [5] 

2. SIMULATIONS 

In this section, we describe the initial power spectra 
and other characteristics of our simulations for the study 
of the evolution of the i/>DM halo MF at high redshifts. 

2.1. Initial Power Spectra 

The suppression of ^DM linear density power spec¬ 
trum relative to CDM can be expressed as 


Pi/jDM{k,z) = T^ BM (k,z)PcDM(k,z), 


( 1 ) 


where P denotes the power spectrum and X^dm is the 
■0DM transfer function (strictly speaking it is the ra¬ 
tio between the transfer functions of ^DM and CDM). 
In general T ^,dm is both redshift- and scale-dependent 
since the balance between gravity and quantum pressure 
introduces a redshift-dependent Jeans scale, kj(z), be¬ 
low which the structures cannot grow. However, for the 
particle masses, redshift range, and halo masses of inter¬ 
est in this work (m 22 ~ 1, z ~ 4— 10, Mh > lx 10 9 -M©), 
Xi/jdm can be approximated as redshift-independent, as 
we demonstrate below. 

The redshift evolution of ipDM density perturbations 
during the matter-dominated epo_ch can be described an¬ 
alytically by (Woo & Chi uehj 2009) 


Pk(k,z) = A(k ) 


3 cos 9 — 9 2 cos 9 + 39 sin 9 
02 ’ 


( 2 ) 


where pu is the spatial Fourier component of the co¬ 
moving density perturbations, A(k) is the normalization 
constant, 9{k,z) = hk 2 \J 1 + z/m^Ho\/D m o, Ho is the 
present Hubble parameter, and fl m o is the present mat¬ 
ter density parameter. Setting 9 2 = 6 gives the Jeans 
scale, 


kj(z) ~ 69.1 to. 


1/2 


Ono h 2 \ 1/A n 
~+iT (1 


-z)- 1 / 4 Mpc -1 , 


(3) 


where h is the dimensionless Hubble constant. For k -C 
kj , we have pk oc (1 4- z) _1 , and thus t/^DM grows like 
CDM; while for k kj the perturbations oscillate as 
Pk oc cos 9. Note that kj(z) increases slowly with time, 
and hence an oscillating mode may become a growing 
mode at lower redshifts, but not vice versa. 

To quantify the difference in growth rate between i)DM 
and CDM density perturbations during a given redshift 
interval, z s tart > z > 2 en d, we define 


m = 


Pk{k"> Zend)/Pk{ ^start) 
Pk{k 0} ^end)/ Pk{kdi 2-start) 


(4) 


with ko <C kj so that the denominator represents the 
growth in CDM. For k <C kj, ^DM grows like CDM, 
and thus £ ~ 1. For k ~ kj, quantum pressure starts to 



Fig. 1.— Growth rate ratio between -0DM and CDM density 
perturbations (Eq. pi) during 30 > z > 4. Vertical dashed lines 
highlight the mass oi the faintest galaxies currently observed by 
HST at z > 4 (M H st) and that expected for JWST assuming a 
limiting absolute magnitude of M\jy = —15 at z = 6 (Mjwst)- 
Here w e have adopted the mass-luminosity relation described in 
Section |4.2.11 Note that even for 77122 — 0.8 the growth rate ratios 
still reach £ ~ 0.99 at = Mhst and £ ~ 0.97 at = Mjwst? 
indicating that the additional suppression of ip DM halos above 
Mjwst during this redshift interval of interest is almost negligible 
(see text for details). 


counter gravity and leads to £ < 1, indicative of ‘addi¬ 
tional’ suppression of ^DM halos during this epoch. 

Figure fTj shows £(fc) for various t/;DM particle masses. 
We take z sta rt = 30 so that the t/jDM density pertur¬ 
bations are still in the linear regime, and z en d = 4 to 
bracket the redshifts of interest in this work. We con¬ 
vert the wavenumber to the halo virial mass via Mh = 
47r(7r/fc) 3 p m /3, where p m is the comoving matter density. 
A relatively large deviation from CDM is found at the 
low-mass end for a smaller particle mass because of the 
corresponding longer Jeans wavelength. However, note 
that the faintest galaxies currently observed by the Hub¬ 
ble Spac e Teles cope (HST) at z > 4 have M h ~ 10 10 M 0 


(see Sec. 4.2.1 for the mass-luminosity relation adopted), 
at which £ ~ 0.99 for to 22 = 0 .8. Even fo r theJames 
Webb Space Telescope (JWST, [Gardner et al. 2006) as¬ 
suming a limiting absolute magnitude of Aluv ~ —15 
at 2 ; ~ 6, we still have £ ~ 0.97 for m 2 2 = 0.8. It 
demonstrates that for the particle masses, redshifts, and 
halo masses of interest when comparing with current and 
forthcoming observations, (i) the growth rate of linear 
density power spectra in CDM and i)DM are similar, 
and (ii) the t/>DM transfer function dm can be well 
approximated as redshift-independent. This is primarily 
due to that the Jeans mass at z = 30 for m 2 2 = 0.8 
is ~ 2.7 x lO 8 Af 0 , well below the observational limits. 
Note, also, that the smallest halos resolved in our simu¬ 
lations have Mh ~ 3 x 10 s M 0 , close to the Jeans mass 
and hence f ~ 0.52 for 771.22 = 0.8. Therefore the halo 
MF at the low-mass end may be, in this sense, slightly 
underestimated in our simulations. 


The UDM transfer function at z = 0 is given by (Hu 
et al.]|000) 


T, 


j/iDM 


cos a; 3 i/is k 

- x = 1.61 to 2 o — 

8 fcj.e, 


1 


(5) 
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Fig. 2.— Linear matter power spectra at z = 30 in CDM and 
^DM. ^DM power spectra are obtained using Equations 0 and 
151, where we have assumed that T- (/ , |, \i can be well approximated 
as redshift-independent during 30 > z > 0 (see Fig. [Tj. Note that 
the smaller the particle mass, the stronger the suppression at the 
high-k end. Arrows indicate the half-mode wavenumbers, fcj /2 (Eq. 
ml ), where the power spectra drop by a factor of four compared to 


DM. 


which is assumed to be redshift-independent relevant 
for the redshifts and wavenumbers of interest. Here 
fcj.eq = 9m22 2 Mpc -1 is the Jeans wavenumber at 
matter-radiation equality. Figure [2] presents the linear 
matter power spectra at 2 = 30, which exhibit sharp 
breaks at k ~ fcj ieq and strong oscillations for k > kj^ eq . 
We can also define the characteristic scale to be the 
‘half-mode’ scale, Aq/ 2 , where T^DMfe/ 2 ) = 1/2- Equa¬ 
tion § then gives 


fci /2 « 5.1 m 2 { 9 Mpc 1 , M 1/2 ~ 3.8 x IO 10 m 22 4/3 Mq, 

( 6 ) 

where M 1 / 2 = 47r(7r/fc 1 / 2 ) 3 p m /3 is the characteristic halo 
mass where a noticeable difference between ^/DM and 
CDM MFs is expected. Note that the strong suppression 
at k ~ fcj, eq shown by Equation ([ 5 ]) is mainl y dete rmined 
during the radiation-dominated epoch (Hu et al.||2000j[, 
and thus cannot be solely explained by EquatiolT 
which is only valid during the matter-dominated epocm 


2.2. Simulation Characteristics 

Bona-fide f/>DM simulations involve solving the 
Sclrrodinger-Poisson equation, which is extremely chal¬ 
lenging owing to its wave nature. The matter wave 
dispersion relation, w oc A~ 2 oc v 2 , where w, A, v are 
the angular frequency, wavelength, and velocity, respec¬ 
tively, indicates that exceptionally high spatial and tem¬ 
poral resolutions are required for resolving the wave func¬ 
tions of high-velocity flows throughout a simulation box 
(SCB14a). Numerically, we find that a comoving spatial 
resolution as high as ~ 1 h~ x kpc is required to prop¬ 
erly resolve a flow with a moderate peculiar velocity of 
~ lOOkms -1 at z ~ 13. Otherwise we find the flow 
velocity can be underestimated, leading to lower mass 
accretion rate and underestimation of MF. In the ex¬ 
treme case, a ^DM simulation in a 30 /z -1 Mpc box with a 
uniform 1 /zA 1 kpc spatial resolution will consume ~ 400 
terabytes of memory, which is impractical in any mod¬ 


ern supercomputer. As a result, even the state-of-the-art 
i/’DM simulations currently can only fully resolve a co¬ 
moving box as small as 1.4/i _1 Mpc to z = 0 (SCB14a). 

In this paper we mainly focus on determining the i/’DM 
halo MF above ~ 1 x 10 9 M 0 at z > 4, for which most 
halos are isolated and thus insensitive to the subtle dif¬ 
ferences between CDM and i/’DM halos. Nor are we in¬ 
terested here in the complex wave nature of the internal 
density profiles of the halos, which we have already estab¬ 
lished in our previous wave-based simulations (SCB14a, 
S14b). As we demonstrated in the previous subsection, 
the growth rates of density perturbations in CDM and 
-i^DM are similar in the context of this work and differ 
mainly in their initial amplitudes. Moreover, S14b veri¬ 
fied that CDM and i/’DM halos have similar virial masses 
during the same collapse process. All these facts indicate 
that, for the purpose of this study, it is appropriate to 
use simulations of collisionless particles with -i/’DM ini¬ 
tial power spectra to approximate real -0DM simulations. 
This is the approach adopted in this work, which is essen¬ 
tially the same as most WDM simulations where initial 
thermal velocity are ignored. Real i/’DM simulations for 
supporting these arguments, solving either wave func- 
tion directly or an alternat ive fluid-like form (e.g., Mocz 
fe Succi|[2015l|Marsh||2015D , are for future work. 


Ail simulations are run from 0 = 100 to 4. Since the 
linear power spectra relevant for this study do not change 
in shape after z = 30, we can directly apply Equation (Dl) 
to obtain the i/’DM power spectra at z = 30 (see Fig. OT. 
To capture the rare non-Gaussian peaks, which are me 
seeds of first galaxies, due to nonlinearity set in as early 
as z ~ 30, we then extrapolate the z = 30 spectra to 
z = 100 for which the amplitude is ~ 3.3 times smaller to 
ensure all perturbations are Gaussian. By doing so, the 
simulations of collisionless particles preserve the shape 
of the spectra from z = 100 to z = 30 but allow for the 
development of rare non-Gaussian peaks. 

We perform simulations with the GADGET-2 N- 
body code (|Spring eI| |200 5|). We adopt the CAMB pack¬ 
age (Lewis et al. 2000[) to generate the CDM trans¬ 
fer function, and construct initial conditions for simu¬ 


lations using the MUSIC code (Hahn & Abel 2011). 


We adopt the cos mological parameters consistent with 
the WMAP9 data (Hinshaw et al. 20131: Q m 0 = 0.284, 
AIao = 0.716, h = 0.696, as = 0.818, and n s = 0.962. 
We choose three different simulation configurations with 
(L,N) = (15 h- 1 Mpc, 512 3 ), (15 h~ x Mpc, 1024 3 ), and 
(30 h -1 Mpc, 1024 3 ), where L is the comoving box size 
and N is the total number of simulation particles. The 
corresponding simulation particle masses are ~ 2.8 x 
10 6 Mq and ~ 3.6 x 10 5 Mq for the lower and higher 
mass-resolution simulations, respectively. For each sim¬ 
ulation configuration, we run four different dark matter 
models: CDM, -i/’DM with m 2 2 = 0.8, 1.6, and 3.2. 

3. MASS FUNCTION 

The main aim of our simulations is to determine the 
halo MF as a function of -i/’DM particle mass. Intuitively, 
a sharp break in the initial power spectrum should trans- 
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Fig. 3.— Density field at z = 4 for t/>DM simulations with 777-22 = 
1.6 in a 15 h _1 Mpc box. Each image displays a projected field for 
a 3 h~ 1 Mpc thick slab with a size of 2.70 x 1.35 h~ 1 Mpc. White 
and red circles show halos more massive than 2 x 10 7 Mq in the 
512 3 and 1024 3 simulations, respectively, where the radii of circles 
equal the halos’ virial radii. The most massive halo has a mass of 
~ lx 10 12 Mq. The upper panel shows both genuine and spurious 
halos, while the lower panel only shows genuine halos. Suspicious 
low-mass halos, which are mostly confined along filaments and have 
no clear counterparts in the 512 3 and 1024 3 runs, are identified as 
spurious, while only massive halos with a good match between low- 
and high-resolution simulations are regarded as genuine. 


dally at low masses (Wang & White 2007 

Angulo et al. 

2013 

Schneider et all] 2013 

)• 


tation due to numerical artifacts (Wang & White 2007), 
and are mostly confined along cosmic filaments (see Fig. 


I 


upper panel). They outnumber genuine halos below a 
laracteristic mass, whic h linearly depends on t he mean 
interparticle separation (|Wang & White||2007|), result¬ 
ing in a prominent upturn m MF at the low-mass end 
(see Fig. [4j open symbols). We define ‘protoh alo’ as 
the in i tial p article positions of an identified halo. |Lovell| 


et al. (2014| showed that the protohalos of genuine and 
spurious halos have distinct features. Genuine protoha¬ 
los are spheroidal and have a good match between low- 
and high-resolution simulations, while spurious protoha¬ 
los have disc-like shapes and their masses and positions 
are sensitive to the simulation resolution, and thus do not 
have clear counterparts in simulations with different res¬ 
olution. To identify and remove these artificial h alos , we 


thus adopt a similar approach suggested by [Lovell et a l. 
(2014) based on the shape of the protohalos and the spa¬ 
tial overlap between low-resolution protohalos and their 
high-resolution counterparts. See Appendix [A] for a more 
detailed description of the algorithms adopted. 

Figure [4] shows the halo MF obtained in our simu¬ 
lations. For comparison, we show both the ‘original’ 
MF (containing both genuine and spurious halos) in 
the 30/i _1 Mpc simulations with 1024 3 particles, and 
the ‘genuine’ MF (with spurious halos removed) in the 



Fig. 4.— Halo mass function (MF) in logarithmic mass bins. 
The open symbols represent the original MF containing both gen¬ 
uine and spurious halos in the 30 h -1 Mpc simulations with 1024 3 
particles, and the filled symbols show the genuine MF with spuri¬ 
ous halos removed in the 15 h -1 Mpc simulations using 512 3 and 
1024 3 particles for estimating the spatial overlap factor. The spu¬ 
rious halos outnumber genuine halos at low masses, resulting in 
an unphysical upturn at the low-mass end of the original MF, es¬ 
pecially for lower 77122 - By contrast, the genuine MF reveals a 
prominent drop at the low-mass end, as anticipated and seen in 
our high-resolution wave-based simulations, reported earlier. The 
shaded regions indicate the uncertainties of genuine MF by varying 
Scut and Ocut by d=20% (see Appendix |A[). Various lines show the 
analytic form, Equation I?}, which fit the simulation results well. 
Arrows mark the minimum t/>DM halo masses proposed by S14b 
for 77122 = 0.8. 


15 ft -1 Mpc simulations using 512 3 and 1024 3 particles 
for estimating the spatial overlap factor. The original 
■0DM MF shows a prominent upturn at the low-mass 
end due to the contamination from spurious halos, es¬ 
pecially for lower TO 22 - By contrast, the genuine i/)DM 
MF features a clear drop at low masses for all redshifts 
and particle masses, apparently different from CDM and 
in agreement with the expectation from a sharp break 
in the ^DM initial power spectra. It is also consistent 
with the minimum ifDM halo mass at z > 1 proposed by 

S14b, M m i n = 3.7 x 10' m .22 ^ (1 + 2) 3 ^ 4 Mq (indicated 
by arrows in Fig. [4] for TO 22 = 0.8). On the other hand, 
the original and genuine CDM MFs are almost indis¬ 
tinguishable, which is no surprise since we assume most 
CDM halos are genuine when calibrating the thresholds 
for removing spurious halos. 

The shaded regions in Figure [4] indicate the uncer¬ 
tainties of genuine i)DM MF by varying S cut and O cut 
by ±20% (see Appendix Aj. It shows that the finding 
of strong suppression of low-mass halos in t/jDM is re¬ 
liable, but the exact slope at the low-mass end is still 
uncertain. In the high-mass end (Mh > 10 11 Mq) the 
original MF is smoother because of a larger simulation 
box. Note, however, that in the intermediate mass range 
(Mh ~ 3 x 10 9 — 1 x 10 11 Mq) the original and genuine 
MFs are reasonably consistent with each other, suggest¬ 
ing that in this mass range (i) most halos are genuine and 
(ii) a 15 F _1 Mpc simulation box is sufficient to obtain an 
accurate MF. These results make the i/jDM MF obtained 
in this work more robust for the purpose of comparing 
with observations. 









































6 


H-Y. Schive et al. 


As shown in Figure [4j the genuine ^DM MF can be 
well fitted by the following analytic form: 


dn 


dM h 


(M h ,z) = 


i/jDM 


dn . /Mh X 11 

®h C DM ( h ’ Z) [ 1+ VMo 

(7) 

where dn/dM h is the halo MF and Mq = 1.6 x 

10 10 tn 22 Mq is the characteristic mass below which 

MF starts to drop noticeably. CDM corresponds to 
77722 —t oo. The facts that Mo has the same particle 
mass dependence as the half-mode mass M-, / 2 (Eq. [iT 
and t/’DM MF drops by a factor of two relative to C 
at Mh ~ Mi /2 reinforce our simulation results. Also 
note that the suppression term, (1 + (Mh/Mg) -1 ’ 1 ) -2 ' 2 , 
is redshift-independent. It is expected since in this work 
the effect of quantum pressure is taken into account only 
for the initial conditions. In detail for t/^DM the suppres¬ 
sion of low-mass halos will be redshift-dependent, but 
the characteristic mass M 0 is still expected to be almost 
redshift-independent since it is mainly determined dur¬ 
ing the radiation-dominated epoch (Hu et ahp OOO). We 
emphasize that the faintest galaxies currently o bserve d 
at x > 4 have M h ~ 10 10 M 0 (see Fig. [B]and Sec. 4.2.1), 
which is close to Mo and hence is insensitive to the un- 
certainties at low masses of MF (M h < 10 9 M 0 ) caused 
by neglecting the dynamical effect of quantum pressure 
and the removal of spurious halos. Equation Q thus 
provides a very convenient comparison between models 
and observations (see next section). 

4. PREDICTIONS VS OBSERVATIONS 

The rest-frame UV LF at high redshifts have become 
increasingly well defined and hence useful for testing a 
range of dark matter models in detail. The physical 
mechanisms assumed to solve the small-scale issues of 
CDM in the Local Group will likely suppress the forma¬ 
tion of faint galaxies at high redshifts as well, so that it 
is important to examine whether too few high-z galaxies 
are created and with a too small Thomson optical depth 
to CMB. This is a well-known iss ue for WDM , usually 
termed as the Catch 22 problem (Maccio et al. 2012). 
In this section we examine the level of consistency in the 
case of (’DM. 

4.1. Cumulative Mass Function 

The cumulative galaxy number density, defined as the 
total number of galaxies per unit comoving volume at 
a given redshift, can be converted into a lower limit of 
m 22 , below which the ^DM MF cannot account for the 
observed counts of galaxies. To relate the UV magni¬ 
tude Myv of a galaxy to its corresponding halo mass, 
we first adopt th e abundance matching technique (Vale 
& Ostriker|2004) which equates the cumulative UV LF, 
4/(<Muv, z), to the cumulative halo MF, n(>A4, z). An 
alternative approach using the conditional LF formalism 
will be discussed later. 

For a given LF, one can apply abundance matching 
to either t/>DM MF with the parti cle masses of interes t 
(Bozek et al. 2015) or CDM MF (Schultz et al. 2014), 
both of which have advantages and disadvantages. The 
former provides a model-independent constraint since 
it simply checks whether the total numbers of t/>DM 



Fig. 5.— Cumulative mass function (MF) at z = 6 — 10. The 
shaded regions indicate the 2a uncertainties. Error bars show the 
observational constraints (2cr at. 2 = 0 — 8 and In at z = 10), which 
match the CDM MFs perfectly simply because they are derived 
from applying abundance matching to the CDM MFs. Note that in 
this approach the faintest galaxies currently observed at z ~ 4 — 10 
all have Mj, ~ 10 ll) Mg. ‘('DM cumulative MFs have finite upper 


and 1.6 are consistent with the observations jBouwens et al. 2007 

McLure et al. 2013 

iSchenker et al. 2013 hinkelstein et al. 20X4 

Oesch et al. 2014) fl 

douwens et al. 
d, t /T _ i nil) 

2015bf, while 77122 = 0.8 has 


details). 

halos at various redshifts are sufficient to account for 
the observed counts, regardless of the underlying mass- 
luminosity ( M-L ) relation. However, the inferred rnass- 
to-light. ratio features a sharp, and probably unphysical, 


drop at the faint end (Bozek et al. 


2015). This approach 


therefore provides a more conservative estimation of 77722 - 
By contrast, the latt er leads to a power -law M-L rela¬ 
tion at the faint end (Schultz et al.l2014), which is more 

fflv 


plausible. However, it fundamentally assumes that CDM 
matches the observed LF perfectly and that ()DM follows 
exactly the same M-L relation as CDM, both of which 
may not be necessarily true. Consequently, any suppres¬ 
sion in the (iDM MF translates directly into a deficit of 
galaxies, resulting in a higher, and likely overestimated, 
lower limit for 77722 - 

Figure [5] shows our cumulative MFs for both CDM and 
7 / 7 DM. CDM MF is constructed from a 30 h~ l Mpc sim¬ 
ulation with 1024 3 particles, and i/jDM MF is obtained 
by integrating Equation 0. The 2cr uncertainties are 
estimated using bootstrap resampling over 125 subvol¬ 
umes, each with a side length of 6/i _1 Mpc. Here the 
halo masses corresponding to various observational data 
points are determined by applying abundance matching 
to the CDM MF, which essentially forces observations to 
be consistent with CDM. In this approach, we find that 
the faintest galaxies currently observed at z ~ 4 — 10 all 
have Mh ~ 1O 1O M 0 . We demonstrate that ^DM with 
77722 = 3.2 and 1.6 are still consistent with the observa¬ 
tions, while 77722 = 0.8 does not have a sufficient number 
of halos with Mh ~ 10 10 M 0 at z = 6 — 8. Using the 
recent LF of B 15b leads to a lower limit of (jDM particle 
mass, 77122 > 1.5 (2cr). By contrast, to infer the results of 
applying abundance matching to the (DM MF, one can 
simply move the observational data points in Figure[5]to- 
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TABLE 1 

Parameters of the Conditional LF Model 


Model 

L 0 

(M A b) 

Mi 

(M 0 ) 

E 

P 

g 

r 

Y 2 

^red 

CDM 

- 20.7 

2.7 x 10 11 

0.16 

1.6 

1.2 

1.9 

1.4 

77122 = 3.2 

- 20.9 

3.1 x 10 11 

0.16 

1.5 

1.2 

1.9 

1.5 

77122 = 1.6 

-21.1 

4.0 x 10 11 

0.16 

1.4 

1.2 

1.9 

1.9 

77122 = 0.8 

-21.7 

7.8 x 10 11 

0.16 

1.2 

1.1 

1.8 

3.1 

B15b a 

-21.9 

1.2 x 10 12 

0.16 

1.2 

1.0 

1.5 



a |Bouwens et al. (2015bJ. 

ward smaller halo masses until, if possible, touching the 
ipDM cumulative MF with a specific particle mass. This 
approach decreases the lower limit to m 2 2 > 0.9 (2a), 
and hence significantly reduces the tension between ob¬ 
servational constraints and smaller ip DM particle masses. 

We emphasize that the two estimations of i/)DM parti¬ 
cle mass given above, namely, 77122 > 1.5 and 777-22 > 0.9, 
are determined using two extreme models of M-L rela¬ 
tion, and therefore likely bracket the uncertainty of the 
lower limit of 77122 - In fact, for t/iDM, it does not make 
much sense to apply the M-L relation predicted by CDM 
from abundance matching. The lower limit 777-22 > 1-5 is 
therefore likely overestimated. In the next subsection we 
provide a more plausible estimation of 77122 based on a 
less model-dependent M-L relation. 



-22 -20 -18 -16 -14 -22 -20 -18 -16 -14 


Wjv [A4b] 

Fig. 6.— Luminosity function (LF) at z = 5 — 10 predicted by 
the conditional LF model. The shaded regions indicate the 2<r 
uncertainties. Error bars show the observed L Fs (2o~ at 2 = 5 — 8 
and la at z = 10) of B15b (open circles) and Oesch et al. (|2014 
open triangles). t/iDM LF shows a drop at the faint end due to the 
suppression of low-mass halos. 77122 = 3.2 and 1.6 are consistent 
with the observations, while 77122 = 0.8 shows a deficit of faint 
galaxies with Muv —18, especially at 2 < 6 (see text for details). 
Vertical dashed lines highlight the limiting absolute magnitudes of 
JWST, which are assumed to be two magnitudes fainter than those 
of HST. 


4.2. Luminosity Function 
4.2.1. Conditional Luminosity Function 

As our preferred method to constrain the t/i DM parti¬ 
cle mass we adopt the conditional LF model (Cooray & 


Milosavljevic 


2005), which has been shown to be able to 


reproduce well the high-z UV LF in the context of CDM 
(Bouwens et al.l 2008 |2015b ). The conditional LF, de¬ 
noted as (p c (L\M\„ z), describes the probability density of 
halos with mass M h hosting galaxies with UV luminosity 
of L. It is modeled by a lognormal distribution, 


(p c (L\M h ,z) = 


1 


V2n\n(10)T,L 


exp 


log[L/Lc(A4,z)] 2 

2E 2 


( 8 ) 

which has a dispersion of ln(10)E and peaks at L C (M\„ z ), 
the M-L relation of the central galaxy. Following B15b, 
we parameterize L c as 


L c (M h ,z) = L 0 


(Mh/MQP f l + z 
1 + (M h /Mi)9 ^ 4.8 


(9) 


where M\ gives the characteristic halo mass. The M-L 
relation asymptotes to L c oc M^ when Mh <C Mi and 
L c oc M^~ q when Mh M\. For a given <p c , the LF can 
then be calculated by 

<t>( L , z ) = <t>c(L\M h ,z)—™-(M h ,z)dM h , (10) 

Jo dM h 

where dn/dM\, is the halo MF. 

Given the above conditional LF formalism, we then 
use chi-square fitting on the observed LF of B15b at z = 
5 — 10 to determine the parameter set (Lq, Mi, E, p, 
q , r). Table 1 shows the best-fit parameters and the 
corresponding reduced chi-square (Xred)- We fix E to 
0.16 both because it is not well constrained owing to the 


substantial uncertainties at the bright end and because 
it does not influence the faint-end slope, which is most 
important for constraining 77122- Note that the faint-end 
M-L relation, L c oc M^, is flatter for smaller 777-22 so 
as to compensate for the stronger suppression of faint 
galaxies. In addition, note that from Equation pi the 
faintest galaxies currently observed at z ~ 4—10 alThave 
Mh ~ 10 10 Mq, consistent with the results obtained by 
applying abundance matching to CDM (see Fig. [ 5 ]). A 
limiting absolute magnitude of Muv ~ —15 at z ~ 6, 
appropriate for JWST, corresponds to Mh ~ 4x 10 9 M 0 . 

Figure [6] shows our predicted galaxy UV LF at z = 5 — 
10 using tire conditional LF formalism described above. 
We use the ST99 MF for CDM, and combine it with the 
ratio given in Equation 0 to get the t/jDM MFs with 
various 777-22- We add 2er variations estimated from the 
MF at Mh ~ 10 10 Mq in our 30/i _1 Mpc simulations to 
capture the uncertainties of the predicted LF around the 
faintest LF bins of B15b. The ^DM LF shows a clear 
decline at the faint end, which is distinctly different from 
the CDM prediction and will be directly testable with 
forthcoming observations such as JWST. Note that this 
feature results from the assumption of a power-law M-L 
relation at the faint end (Eq. [ 9 ]), and thus cannot be 
captured by the usual abundance matching. The cases 
77722 = 3.2 and 1.6 are found to be consistent with the 
current observations, while 77122 = 0.8 shows an apparent 
deficit of faint galaxies with Muv A) —18, especially at 
z < 6. This result is consistent with the analysis of 
cumulative MF (see Figure [5]). Us ing the faintest LF 
bins of B15b (open circles in Fig. |6| leads to 77122 > 
1.2 (2tr). Note that thi s constraint lies in-between those 
obtained in Section 4.1 (77722 > 0.9 and 1.5), as expected, 
since the faint-end M-L relation adopted here can be 
regarded as a compromise between the two extreme cases 
using abundance matching to CDM and to t/jDM MFs, 
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Fig. 7.— Luminosity function (LF) at z = 2 — 4 predicted by 
the conditional LF model. The shaded regions indicate the 20- 
uncertainties. The data with 2 o error bars are the LFs determined 
from B15b (open circles) and |Parsa et al. (2015] open squares). 
The observed faint-end slope m this redshift range is consistent 
with the i/’DM model with 17122 ~ 1.6 (see text). Notice that CDM 
overestimates the number density of faint galaxies, especially at 
2 = 3. 


respectively. 

It is also interesting to extend the comparison of LF 
to z < 5 (Fig. [7]). To be self-consistent, we adopt the 
same pa rameters o f the con ditional LF model shown in 
Table 1. Parsa et al. (2015) recently found a much shal- 


from the steep slope reported previously (Reddy & Stei- 

del 

O 

O 

CM 

Alavi et al. 2014 

Bouwens et al. 2015b). t/jDIVI 


this shallower slope, while CDM overestimates the num¬ 
ber density of faint galaxies, especially at z = 3. It is 
thus very important for future research to understand 
this apparent discrepancy betwe en the faint-end slopes 
determined by Parsa et al. (2015) and the previous stud¬ 
ies at z < 5. 

It should be emphasized that the conditional LF model 
provides a more reasonable constraint on the i/jDM parti¬ 
cle mass. At first glance, it may seem that this approach 
has too many free parameters to provide an appropri¬ 
ate estimation. However, actually the only relevant pa¬ 
rameter for constraining 77122 is p , the faint-end slope 
of the M-L relation, as argued below. Since the faint- 
end slope of the LF is insensitive to S, we have L ~ 
L c oc MThis leads to </>(Muv) = 0.4ln(10)L^(L) = 
0.41n(10)p _1 dn/dlnMh, where dn/dlnM h is the halo MF 
in logarithmic mass bins. Therefore, for a given p, a 
maximum observed </>(Muv) can be directly converted to 
a minimum required peak dn/ dlnMh, which then turns 
into a lower limit of 77122 - Moreover, if dn/dlnM h oc MJL 
which is appropriate when M h Mq in Equation (Ffl 


(i.e., when -0DM is still close to CDM), then <f>{L) oc 
L'n/p~ l . Since p can be estimated by the ST99 MF, we 
can determine p from a given <j>(L). Accordingly, p itself 
is not unconstrained. 

As an example, consider the t/iDM model with 77722 = 
1.6 at z = 6. The observed LF has oc L~ 21 at 
Muv ~ —19. This luminosity corresponds to a halo mass 
of Mv, ~ 7 x 10 10 Mq, at which dn/dlnM h oc M/f 1 * ' 5 * 7 * * (see 
Fig. [4]). We thus have p = —1.5/(—2.1 + 1) ~ 1.4, con¬ 
sistent with Table 1. The MF has a peak of dn/dlnM h ~ 
(1.0 ± 0.3) x 10 -2 Mpc -3 around M^ = 10 10 Mq, which 
can then be converted to a maximum LF of </>(Muv) = 
(6.8 ± 2.0) x 10~ 3 Mpc _3 . This estimation is in excel¬ 
lent agreement with Figure [6l even though the only as¬ 
sumption made here is that the M-L relation is a power 
law at the faint end. Therefore, we conclude that the 
conditional LF model provides a more plausible and less 
model-dependent estimation of the t/’DM particle mass 
than with abundance matching. 


4.2.2. Truncated Schechter Function 

It is also useful and convenient to parameterize the pre¬ 
dicted LF of -0DM by a formula similar to the Schechter 
function. We adopt 


m = 


L* 



where 


F(L) 


1 + 



( 11 ) 


( 12 ) 


represents the suppression of faint galaxies in ^>DM (F = 
1 for CDM). L^ is the characteristic luminosity of the 
suppression, below which cj>(L) asymptotes to L“ +/3 . To 
describe the time evolution of LF, we follow the literature 


3ouwens et al. 2012 Kuhlen fe Faucher-Giguere 2012 


Schultz et al. 2014 Bouwcns et al. 1 2015a) and assume 


that the parameters in Equations (TT) and (121 depend 
linearly on redshift. Applying chi-square fitting to the 
conditional LF model (Fig. [6| then leads to 


M uv ,* = -20.90 - 0.004(2 - 6) 

0* = O.52 x io~°- 28 ( z-6 ) -3 Mpc -3 

a = -1.78 -0.06(2 -6) 

Muv,* = -17.44 + 5.19 log(m 22 /0.8) - 2.71 log((l + z)/7) 
/3 = 1.69+ 0.03(2-6) 

7 = —1-10, (13) 


where we have assumed Muv = 

—2.5 log(L/ergs -1 Hz -1 ) + 51.6. 

Figure [8| shows the LF obtained by the truncated 
Schechter function described above. It fits well with the 
LF predicted by the conditional LF model at 2 = 4 — 8, 
while for the bright end at 2 = 10 it slightly outnumbers 
the observed galaxies (see B15b and references therein) 
and is marginally consistent with the conditional LF 
model. This subtle discrepancy mainly results from the 
assu mpt ion of linear dependence on redshift in Equa¬ 
tion (13), and it is unclear whether it indicates a faster 
evolution at 2 > 8 given the substantial uncertainties in 
the LF at 2 = 10. This possible acceleration in evolution 
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Fig. 8 .— Luminosity function (LF) at 2 : = 4 — 10 obtained by 
a singl e analytic formula similar to the Schechter function (Eos. 
|11|13|; central lines). The shaded regions are the same as Fig. |6] 
snowing the LF predicted by the conditional LF model within za. 
E rror bars show the o bserved LFs (2 a at z = 4—8 and lcr at 2 : = 10) 
of | Par sa ct al. (2015| open squares), B15b (open circles), and|Oesch| 
[etahl ( [2Ul4| ~ open triangles). The analytic formula well reproduces 
the conditional LF results at 2 : = 4 — 8, while at 2 : = 10 it slightly 
outnumbers the observed galaxies and is marginally consistent with 
the conditional LF model. 


has been successfully modeled recently in the context of 
abundance matching for CDM by ad ditionally incorpo¬ 
rating early stellar evolution ( Mason et al.}2 015). It may 
be interesting to apply this approach in the context of 
t/iDM too for extending predictions to z > 10 with some 
security. 


4.3. Magnification Bias for the Hubble Frontier Fields 

We have seen above that the quantum pressure inher¬ 
ent to ipDM leads to a suppression of low-mass galax¬ 
ies and hence our predictions for the LF have largest 
contrast with CDM at low luminosities. Here we exam¬ 
ine the benefits of gravitational magnification in the new 
HFF data, where statistical samples of mult iply lensed 
galaxies are magnified by typically ~ 10 (Lam et al. 
20141, reaching two magnitudes or more further down 
the LF at high redshifts. It corresponds to an intrinsic 
UV luminosity of Mjjy ~ —15, where we predict sizeable 
difference between CDM and ^DM for the interesting 
range of m 22 needed to provide the kpc-scale dark cores 
of local dSph galaxies. 

Gravitational lensing induces a bias in the number den¬ 
sity of sources detected above a flux limit, which in the 
case of lensing clusters is well establish ed (Broadhurst 


et al. 1995 Urnetsu & Broadhurst 


the detection of the highest re dshift and lowest lumi- 


2008) and has led to 


nosity galaxies currently known ( Zheng et al.||2012 Coe 


et al. 2013 Zitrin et al. 20141. The number density of 
such high-redshift galaxies is modified in the following 
way: 

Vlensed(^L) = (l/^l) V un i ense d(>Z//p,), (14) 


where ^ is the magnification factor and N(>L) is the 
galaxy number density above a flux limit corresponding 
to L. It shows the competition between the enhanced 
number density due to the lower, magnified limiting lu¬ 
minosity and the diminished source plane area, which is 


Fig. 9.— Magnification bias predicted by 1 / 1 DM and CDM at 
z = 10. The numbers on the top of each panel indicate the adopted 
limiting luminosities. For CDM the magnification bias continually 
rises owing to the steep faint-end slope of the LF, whereas for 
t/)DM the bias is generally lower than one because of the strong 
suppression of low-mass halos. The lower the limiting luminosity, 
the greater the contrast between xj )DM and CDM. 


smaller than the observed area by the same magnifica¬ 
tion factor. The magnification bias can thus be defined as 
Mensed(>£)/-^uniensed(>£)j which equals one when there 
is no bias. 

Figure [9] shows the difference between the magnifica¬ 
tion bias for CDM and -0DM at z = 10, which we predict 
to be a strong function of limiting luminosity because of 
the difference in sign in the faint-end slope of the LFs be¬ 
tween these two models. CDM is continually rising and 
hence the magnification bias is greater than one, thereby 
enhancing the number of faint galaxies detected at high 
redshifts, whereas for i/iDM the turnover in the LF leads 
to many fewer galaxies magnified above the flux limit. 
The difference we predict at M uv = —14.5 for example 
is a factor of ~ 10 for /r ~ 10 at z ~ 10 between CDM and 
■0DM with to 22 ~ 1.2. Going beyond this with JWST 
should probe another two magnitudes d eeper with the as¬ 
sistance of lensing (Mason et al. 2015) - most efficiently 
by employing the same deep lenses as the HFF for which 
the magnificatio n maps have been wid ely studied flRod- 
|ney et ah1|2015|) a nd best understood (|Lam et al.| 2UI4| 
|Diego et al.||2015 1. 


4.4. Reionization 

Based on the predicted -0DM LF, we can calculate the 
reionization history using the standard approach that 
has been adopted for various dark matter models (e.g., 
Kuhlen fe Faucher-Giguere| |2012[ Schultz et al.| 2014[ 
Bozek et al. 2015|). The time evolution of the volume 
filling fraction of ionized hydrogen, Qmi(z), is governed 

by 

dQ HII fhon Q HII /, 

= -T—’ ( 15 ) 

dt 77 »h tree 

where ?Th is the mean comoving hydrogen number den¬ 
sity. We take Qhii =_0 at z = 25. The volume-averaged 
recombination time, t lec (z), can be determined by 


tr, 


0,3 (%iy 


l + z 


-3 


Gyr, (16) 
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where Chii = (Mn)/( n Hii) 2 is the volume-averaged 
clumping factor. Here we have assumed an intergalac- 
tic medium temperature of 2 x 10 4 K and a primordial 
hydrogen mass fraction of 0.76. 

The comoving ionizing emissivity, 77i on (z), defined as 
the number of ionizing photons produced per unit time 
per unit comoving volume, can be estimated from the 
galaxy UV LF: 


TABLE 2 

Reionization parameters 


Model 

Chii 

Cion J esc 

Ml im a 

MIN 

4.0 

0.2 

-13 

FID 

3.0 

0.6 

-13 

MAX 

2.0 

1.0 

-13 


a -^lim is allowed to vary when computing r e . 


2 x 10 25 
erg Hz -1 


Cion fei 


dL(t>(L)L , (17) 


il i in 


where Cion represents the efficiency of converting galaxy 
UV luminosity to ionizing photon luminosity, and / esc 
is the escape fraction. Strictly speaking, since the ob¬ 
served rest-frame UV luminosity (at ~ 1500 A) will also 
be extinguished by dust, f csc in Eq uation ( |17| ) is the rel¬ 
ative escape fraction (Steidel et al]'2001) denned as the 
absolute escape fraction (the traction of ‘ionizing pho¬ 
tons’ that escapes the galaxy without being absorbed by 
dust and neutral hydrogen) divided by the fraction of 
‘UV photons’ that escapes. This relative escape fraction 
can be significantly higher than the absolute escape frac¬ 
tion because of the efficient dust extinction at ~ 1500 A. 
Mum = — 2.51og(Lii m /ergs -1 Hz -1 ) +51.6 is the limit¬ 
ing UV magnitude, below which the galaxy formation is 
assumed to be inefficient. Note that hi on is sensitive to 
the faint-encl slope of LF that differentiates various dark 
matter models. 

For a given Qym(z), the Thomson optical depth to 
CMB can be calculated via 


r e = CCTt n a 


’ n + z ) 2 

dz [ -^-Q mi (z)(l + v(z)Y/4X), 
H(z) 


(18) 

where cr x is the Thomson cross-section, H(z) is the Hub¬ 
ble parameter, X ~ 0.76 and Y = 1 — X are the primor¬ 
dial mass fraction of hydrogen and helium, respectively, 
and we take i](z > 4) = 1 when helium is only singly 
ionized and rj(z < 4) = 2 when helium is doubly ionized 
by quasars. 


There are three free parameters in Equations (15)- 
(18), namely, Chii, XI , Con ./Ac (Cion and / esc are 
fully degenerate), which, for simplicity, are assumed to 
be spatially uniform and redshift-independent in this 
work. The typical parameter ranges adopted in the lit¬ 
erature are Chii = 2 ~ 5, Mn m = —17 ~ —10, Con = 
0.5 ~ 2.0, and / esc = 0.1 ~ 0.5 (Bouwens et al.| 2012 


Kuhlen & Fauchcr-Gi guere| |2012[ |Schultz et al.| 2014 

Bozek et al.|2015[|Bouwens et al.|2015a[ Robertson et al. 


2015|). To bracket the uncertainties in these parameters 

we consider three different parameter sets: a minimum 
reionization model (MIN) with (Chii, Con/esc, Miim) = 
(4.0, 0.2,—13), a fiducial reionization model (FID) with 
(Chii, Con/esc, Mim) = (3.0, 0.6,-13), and a maximum 
reionization model (MAX) with (Chii, Cion/esc, Miim) = 
(2.0,1.0, —13). The adopted values are also shown in Ta¬ 
ble 2. Note that since i/)DM should be insensitive to My lm 
(because of the strong suppression of small halos), we fix 
Mim = ~13 unless otherwise specified. 

Figure [TO] shows the time evolution of Qhii for various 
models. For a given reionization parameters, i(DM pre¬ 
dicts a faster increase of ionized volume at later times. 
Recent observations indicate that reionization undergoes 



Fig. 10.— Volume filling fraction of HII as a function of redshift. 
The central lines correspond to the FID model and the shaded 
regions are bounded by the MIN and MAX models. The limiting 
UV magnitude is fixed to M\[ m = —13. Open ci rcles and arrows 


Schroeder et al. 2013 |Schenker et al. 2014; McGreer et al. 2U15|. 

a rapid evolution at 6 < * < 8 and 

is completed at z ~ 6 

1 Fan et al. 2006 Schroeder et al. 

2013 

Schenker et al. 


quired ionizing photon production efficiency, Con/esc, is 
~ 0.6 for (iDM with 77122 = 0.8, about two times higher 
than CDM (Con/esc ~ 0-3). For CDM this result is in¬ 
sensitive to the values of both M\\ m and Chii probed in 
this work. 

Figure [TTl shows the Thomson optical depth as a func- 


llUll U1 i Klli m VctllUUb IllUUtUO. IVlUoL It 

2015 results give r e , = 0.066 ± 0.016 (1<t, 

Planck Col- 

lab oration et al. 

2015), which is marginal. 

y consistent 


with the estimation from Quu(z). Note that this estima¬ 
tion is largely independent of both Chii and My lm when 
Mum > —15, since galaxies fainter than this magnitude 
are highly suppressed at z > 6 even for 77122 = 3.2. Given 
the considerable uncertainties in the reionization model 
at high redshifts, it is therefore clear that neither Qmi{z) 
nor r e provides a stringent constraint on i/DM even for 
77722 ~ 0.8. In comparison, CDM with an ionizing pho¬ 
ton production efficiency as low as Cion/esc ~ 0.2 can be 
consistent with Planck 2015 as suming My, m > —13, in 


agre ement with the f inding s of Bouwens et al. (2015a) 


and Ro bertson et al . (2015). 

Bozek et al. l2U15p used a similar approach to estimate 


the ipDM particle mass. They reported that 77122 = 10 
is disfavored by the observed value of r e at 3a, and con¬ 
sequently a particle mass as high as 77122 = 10 is re¬ 
quired, significantly higher than our estimation. This 
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Mim [Mab] 


Fig. 11.— Thomson optical depth versus limiting UV magnitude. 
The central lines correspond to the FID model and the shaded 
regions are bounded by the MIN and MAX models, where Mi; m is 
allowed to vary. Th e cross-h atched region sh ows the Planck 2015 
lcr confidence limit (I'lanck Collaboration et al.|2015]|. 


inconsistency with our result mainly arises from a higher 
value of r e they adopted from the p revious Planck 2013 


results, t b = 0.090 ± 0.013 (l cr, Planck Collaboration 
et al. 2014 Spergel et al.||2015| . Also note that the LF 


we predict based on the ‘conditional LF model’ declines 
at the faint end (see Fig. [6|, while the LF adopted by 


Bozek et al. (2015) based on the ‘abundance matching’ 
approach does not have this natural feature. Such rela¬ 
tive deficit of faint galaxies in our LF model would only 
delay the reionization process and reduce the Thomson 
optical depth, hence increasing the discrepancy between 
observations and the t/jDM model with a small parti¬ 
cle mass. However, here we demonstrate that 777-22 ~ 0.8 
can still be consistent with the latest Planck observations 
provided that the ionizing photon production efficiency 
is sufficiently high. 


5. DISCUSSION AND CONCLUSION 

In this paper, we have constructed cosmological simu¬ 
lations designed to study the dark matter halo MF in the 
wave dark matter (0DM) scenario. Here the uncertainty 
principle counters gravity below a Jeans scale, which is 
determined by the only free parameter in this model, 
77722 , the dark matter particle mass. The smaller the 
particle mass, the larger the Jeans scale, and hence the 
stronger the suppression of low-mass halos. For this rea¬ 
son, we focus on determining a lower limit of 77722 based 
on the observed UV LF at z ~ 4—10 and the reionization 
history. 

The major findings in this study can be briefly sum¬ 
marized as follows: 


• i/'DM halo MF has a prominent drop below ~ 
10 10 Mq, which can be well fitted by Equation ([ 7 ]). 

• i/)DM predicts a clear drop in the galaxy LF around 
IlLuv ^ —16 at z > 4 based on a conditional 
LF model, which can b e fitte d by a truncated 
Schechter function (Eqs. EH)- 

• The newly established LF at z ~ 4 — 10 constrains 
the t/iDM particle mass to be 77122 > 1.2 (2cr). 


• For galaxies magnified >10x in the Hubble Fron¬ 
tier Fields, ipDM predicts an order of magnitude 
fewer detections than CDM at z > 10 down to an 
intrinsic UV luminosity of M\jy ~ —15. 


• t/jDM with TO 22 > 0.74 can satisfy the Thomson 
optical depth reported by the latest Planck obser¬ 
vations, on the assumption of a reasonable ionizing 
photon production rate. 


In the following we give a more thorough discussion of 
this work. We first argue that, for studying the t/’DM MF 
with the particle masses, redshift range, and halo masses 
of interest here (77122 ~ 1, z ~ 4 — 10, M j, > lx 10 9 Mq), 
it is reasonable to approximate both the t/iDM trans¬ 
fer function as redshift-independent and the evolution 
of quantum fluid by collisionless particles. The major 
drawback of these approximations is the inability to cap¬ 
ture the difference in the substructures between CDM 
and f)DM halos, where t/>DM halos have prominent soli- 
tonic cores in the centers surrounded by fine-scale, large- 
amplitude cellular interference (SCB14a, S14b). How¬ 
ever, this shortcoming is irrelevant when one is only con¬ 
cerned with the halo masses, as in this work. This is 
especially true because most halos above ~ 10 9 Mq have 
yet to merge gravitationally with each other at z > 4. 

Simulations of collisionless particles with a cut-off in 
the initial power spectrum suffer from the well-known 
side effect of inducing spurious halos due to numerical 
artefacts, which must be accounted for when determin¬ 
ing an accurate MF at low masses. To identify and re¬ 
move these artificial halos, we a dopt a similar approach 
suggested by Lovell et al. (2014) based on the shape of 
the progenitors of halos and the spatial overlap between 
low-resolution halos and their high-resolution counter¬ 
parts (see Appendix 0. The resulting MF cleaned in 
this way then shows a clear decline as expected below 
the Jeans mass and can be well fitted by Equation ([7|. 
Most importantly, this reinforces the MF we derive above 
> 3 x 10 9 Mq, which are most relevant for comparing 
with observations. 

Comparing the halo MF with the observed UV LF re¬ 
quires the knowledge of M-L relation. For this purpose, 
we first apply the abundance matching technique to ei¬ 
ther CDM or t/jDM MFs. In both cases, i/)DM with 
77722 = 0.8 shows a deficit of low-mass galaxies with 
Mh ~ 10 10 Mq at z = 6 — 8 (see Fig. [ 5 ]), leading to lower 
limits of particle mass, 77722 > 1.5 (abundance match¬ 
ing to CDM) and 777-22 > 0.9 (abundance matching to 
^DM). We also explore the conditional LF model as an 
alternative approach, which yields 77722 > 1.2. The key 
assumption here is a power-law M-L relation at the faint 
end. We argue that this approach provides a more rea¬ 
sonable and model-independent estimation of the t/jDM 
particle mass. In addition, we predict that the high-z LF 
should turn over slowly around Muv > —16 at 2 > 4, 
distinctly different from CDM. This predicted feature lies 
just beyond the detected luminosity range of the current 
LFs at z > 4, but will be directly testable with forth¬ 
coming observations such as JWST and also with highly 
magnified galaxies in the HFF data. 

Note that the M-L relation for low-mass halos may 
be subject to large uncertainties. Strigari et al. (2008) 
showed that the Milky Way dwarf satellites, which have 
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Mh < 10 10 Mg, share a common mass scale but have 
luminosity differences over four orders of magnitude. 
Therefore, a more complicated M-L relation might be 
expected for halos below 10 10 M< q, at least at lower 
redshifts. Even in CDM simulations, O’Shea et al. (20151 
found a relatively flat high-z LF at the faint end com¬ 
pared with the Schechter function fits to observations 
(B15b). Ideally, t/>DM simulations with the addition of 
baryonic physics may be very helpful in further differen¬ 
tiating the high-z LFs predicted by CDM and ^DM. 

The Thomson scattering optical depth to CMB pro¬ 
vides another constraint for the i/)DM particles mass. In 
general, i/)DM predicts a faster increase of ionized vol¬ 
ume at later times due to the suppression of early galaxy 
formation. We demonstrate that (/'DM with 77122 > 0.74 
can satisfy the Planck 2015 results and have reioniza¬ 
tion completed at 2 > 6, on the assumption that the 
ionizing photon production rate is sufficiently efficient 
(about three times higher than that required for CDM). 
This result is largely independent of the limiting luminos¬ 
ity and the faint-encl slope adopted since galaxies fainter 
than Muv ~ —15 are highly suppressed at z > 6 even 
for 77722 = 3.2. On the other hand, this constraint is 
somewhat undermined because of the current large un¬ 
certainties associated with the escape fraction and the 
efficiency of converting galaxy UV luminosity to ionizing 
photon luminosity. Note that for CDM the reionization 
calculation is much more uncertain in fact than for 0DM, 
as it is dominated by the choice of limiting luminosity as¬ 
sumed for the relatively steeply rising LF (a > 1.0) as 
otherwise the integrated luminosity diverges. 

The MF below ~ 10 9 Mq determined by our simu¬ 
lations has a larger relative uncertainty in our model. 
Firstly, approximating the )/)DM transfer function as 
redshift-independent can somewhat underestimate the 
matter power spectra at higher redshifts (see Fig. fTh , 
which in turn leads to a underestimation of MF at low 
masses. Secondly, we approximate the evolution of quan¬ 
tum fluid by collisionless particles. It may allow for the 
formation of a small number of halos with masses well 
below the Jeans mass, which would otherwise be sup¬ 
pressed further if the dynamical effect of quantum pres¬ 
sure is taken into account. In this sense, the MF in a 
bona-fide wave-based ()DM simulation may have an even 
stronger break at the low-mass end than Figure [4| We 
may hope to check on this in the future by running the 
full adaptive-mesh-refinement wave simulations that we 
have previously described (SCB14a, S14b) on a substan¬ 
tially more powerful platform. Thirdly, there are also 
uncertainties associated with the removal of spurious ha¬ 
los below ~ 10 9 Mq (see the shaded areas in Fig. [ 4 J. 
However, it should be emphasized that none of these un¬ 
certainties are relevant for the purpose of this study at 
the level of precision that is currently afforded by the 
data, as the depths of the current Hubble and forthcom¬ 
ing JWST do not extend to such low-mass halos. 

Though not entirely consistent, the i)DM particle mass 
estimated in this work, 777-22 >1.2 (2a), is surprisingly 
close to the values determined from local dwarf galaxies. 
SCB14a established with the first wave-based 7/1DM sim¬ 
ulations a distinct solitonic core in the center of every 
i/)DM halo. We have previously obtained to 22 = 0.8±0.2 
(ler) by fitting the spatial distribution of the interme¬ 
diate metallicity stellar population in the Fornax dSph 


galaxy to the soliton mass profile, under the assump¬ 
tion of a con stant projected velocity dispersion. Marsh 
& Pop (2015) have also determined a similar constraint, 
77122 < 1.1 (2 a), by fitting the mass profiles of Fornax and 
Sculptor dSph galaxies to the soliton mass profile using 
an empirical relation between the half-light radius and 
velocity dispersion. Note that for to 22 = 1.2, we predict 
that a halo with Mh = 2 x 10° Mq still has a core as large 
as ~ 1.1 kpc (S14b), consistent with many estimates of 
the large cores found in dSp h galaxies (e.g., Salucci et al. 
2012| |Amorisco et al.||2013 ). A more thorough compari¬ 
son between the stellar phase-space distribution in dSph 
galaxies and the (jDM halo mass profile using a full Jeans 
analysis will be extremely important to further clarify 
how coincident are the particle masses determined by 
these various approaches and the role of baryonic feed¬ 
back in this context (S-R. Chen et al., in preparation). 
The inherent density granularity of the t/jDM halo may 
also be examined through internal dynamical effects and 
by lensing flux anomalies on sub-kpc scales, which pro¬ 
vide other key independent observational tests for distin¬ 
guishing i/jDM from WDM and CDM. 

The finding that a similar particle mass in (/)DM 
can b oth solve t he cusp-core problem in dwarf galax¬ 
ies (Moore 1994) and satisfy the high-z LF and reion¬ 
ization observations is very encouraging for this model. 
It demonstrates a great advantage of t/’DM over WDM, 
for which the light particle mass required for creating 
kpc-scale cores in dSph galaxies prevents the formation 
of the host dwarf galaxies in the first p l ace and over ly 
suppresses high -z galaxies (Maccio et al. 2012 Schnei¬ 
der et al. 2014). The key reason for this striking ditler- 
ence is that the relation between core radius and power 
spectrum truncation differs because of the very differ¬ 
ent physics underlying t/jDM and WDM, with the uncer¬ 
tainty principle responsible in the case of t/iDM and free 
streaming in the case of WDM. The particle masses in 
the two models can always be chosen so that they have 
similar truncated matter power spectra, but the corre¬ 
sponding core radius in ^>DM can then be several times 
larger than that in WDM. This is why t/jDM does not 
suffer from the Catch 22 problem affecting WDM. 
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APPENDIX 

A. REMOVAL OF SPURIOUS HALOS 

Particle simulations with an initial power spectrum 
cut-off suffer from the formation of spurious halos, es¬ 
pecially at low masses. It is therefore necessary to adopt 
a robust algorithm to identify genuine halos in the sim¬ 
ulations, which we describe below. 


We first use the AMIGA Halo Finder (AHF, Knoll- 
rnann & Knebe 2009) to identify all halos (both genuine 
and spurious halos), and find the corresponding proto¬ 
halos by tracing the constituting particles back to the 
initial redshift. We then define the protohalo sphericity 
as 


S = 


I\ ± 1-2 ~ I 3 


-h 


h 


(Al) 
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Fig. Al.— Sphericity and spatial overlap factor as a function 
of halo mass in the 15 /i —1 Mpc simulations. The central solid and 
dotted lines show the mean values, and the shaded areas contain 
95% of halos. Low-mass -0DM halos have lower sphericity and 
overlap factor compared with CDM halos, especially for Mh < 
10 9 Mq, indicative of a severe contamination from spurious halos 
in this mass range. 


where I\ < I 2 < I3 are the principle moments of inertia 
of a protohalo. This is equivalent to having S = c/a , 
where a and c are the maximum and minimum semi-axis 
lengths of an ellipsoid of uniform density with the same 
princ iple moments of i nertia as the protohalo (Lovell 
et al.{2014| . Figure I aTI ( upper panel) shows the spheric¬ 
ity as a function of halo mass for both CDM and >/)DM 
with 17122 = 1-6 in the 15 h~ l Mpc simulations with 512 3 
particles. Clearly, CDM halos have higher S. The av¬ 
erage S of ('DM halos lies below the 2a lower limit of 
CDM halos for M h < 10 9 Mq, indicative of a severe con¬ 
tamination from spurious 0DM halos in this mass range. 

To quantify the matching accuracy between two pro¬ 
tohalos in low- and high-resolution simulations, respec¬ 
tively, we introduce a spatial overlap factor, 


O = 


Vow n Vhigh 
Vow U Vhigh ’ 


(A2) 


where V is the volume of a protohalo estimated by 
depositing particles onto grids using the cloud-in-cell 
scheme with a cloud size equal to the mean particle sepa¬ 
ration. Vow H Viigh and Vow U Vhigh represent the intersec¬ 
tion and union of two volumes, respectively. Accordingly, 
O = 1 if the spatial distribution of two protohalos are 
completely overlapped, while O = 0 if there is no overlap 
at all. Note that our definition of sp atial overlap factor 
differs fro m th at of Lovell et al. (2014 see Eq. 6 therein). 
Equation (A21 does not require solving the gravitational 
potential and thus is more computationally efficient and 
more sensitive to th e ac tual spatial overlap between two 
protohalos. Figure Al (lower panel) shows the spatial 
overlap factor as a function of halo mass. The average O 
of (DM halos shows an apparent drop at Mh < 10 9 Mq, 
again suggesting that spurious halos start to outnumber 
genuine halos. 

We consider most CDM halos as genuine since its MF 
is consistent with the analytic prediction of ST99 above 
~ 3 x 10 8 Mq (corresponding to ~ 100 particles in our 
lower resolution simulations). We thus define the 2a 
lower limit of S and O at Mh = 3x 10 s Mq in CDM as the 
minimum thresholds, S cu t and O cut , below which halos 
are regarded as spurious. S cu t and O cut slightly increase 
with redshift, and we take 5 cu t = 0.31 and Ocut = 0.34 
by averaging over z = 4 — 10. Not e that we do not us e 
the mass cut criterion adopted by Lovell et al. (2014), 
which relies on extrapolating the low-resolution results 
to high-resolution runs. 

In short, for each low-resolution ha lo identified by 
AHF, we calculate S by Equation (Al) and det ermine 
the highest value of O by applying Equation 1A2h for all 
nearby high-resolution halos. Then, we use CDM halos 
to determine the minimum thresholds, S cu t and 0 C ut> and 
only low-resolution halos with S > S cu t and O > 0 C ut 
are marked as genuine. Figure pM (lower panel) shows the 
spatial distribution of genuine Tialos. Clearly, suspicious 
low-mass halos, most visible along filaments, are removed 
with the above procedure, while only massive halos with 
a good match between low- and high-resolution simula¬ 
tions are retained. 
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